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Abstract: In the first part of this paper, we extend the ci-dimensional unitarity cut method of hep- 



ph/0609191 to cases with massive propagators. We present formulas for integral reduction with which one 
can obtain coefficients of all pentagon, box, triangle and massive bubble integrals. In the second part of this 
paper, we present a detailed study of the phase space integration for unitarity cuts. We carry out spinor 
integration in generality and give algebraic expressions for coefficients, intended for automated evaluation. 
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1. Introduction 

With the approach of the Large Hadron Collider experiments, accurate descriptions of particle physics 
will require knowledge of one-loop cross sections. Computational complexity increases dramatically with 
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the number of legs, even at the amplitude level. It is desirable to find a simple and fast algorithm for these 
computations. There has been notable progress in the last couple of years [|l], ||, ||], [j| |(| 0, ||, ||, [To|, 



The unitarity method introduced in 12] seeks to compute amplitudes by applying a unitarity cut 
to an amplitude on one hand, and its expansion in a basis of master integrals on the other [13]. With 
knowledge of the basis and the general structure of the coefficients in the expansion, the coefficients can 
be constrained. 



The holomorphic anomaly [14] reduces the problem of phase space integration to one of algebraic 
manipulation, namely evaluating residues of a complex function. By applying this operation within the 
unitarity method, coefficients can be extracted systematically. In this manner, a method was introduced 
to evaluate any finite four-dimensional unitarity cut and systematically derive compact expressions for the 
coefficients || . 

When working with four-dimensional cuts, one loses information of possible rational terms that are 
cut-free in four dimensions. However, unitarity methods can find rational terms as well, if we carry out the 
integral in d = (4 — 2e) dimensions to higher orders in e [15]. This program was developed in |l6|, 17, 15, 

Recently, a general d-dimensional unitarity cut method was developed for one- loop amplitudes [19, 20]. 
In this method, coefficients are extracted by first separating and performing a four-dimensional integral 
by a technique of choice, and then identifying the integral over the remaining d — 4 dimensions with a 
coefficient using recursive dimensional shift identities. In principle one can work out the coefficients to all 
orders in e. In fact, for cases with only massless propagators, this method is a complete alternative to 
Passarino-Veltman reduction. 

In the first part of this paper (Sections 2 and 3) we generalize the work of (l^, |2(| to cases where 
propagators have nonzero, non-uniform masses. We present formulas for integral reduction from which 
one can obtain coefficients of all scalar pentagon, box, triangle and massive bubble integrals. Pentagons 
in arbitrary dimensions are truly independent master integrals, though in four dimensions they can be 
expressed in terms of boxes up to 0(e). Amplitudes allowing massive propagators will most generally have 
tadpole and massless-bubble master integrals as well. Unitarity methods cannot access this information 
directly, but other physical considerations pin it down. 

In the second part (Section 4), we further simplify the technique of spinor integration over phase 
space by carrying out some intermediate steps in generality. These steps include expressing the integrand 
as a total derivative and a possible Feynman-parameter integral. We give algebraic expressions for the 
coefficients of the cut-containing basis integrals. Our formulas should be suitable for programming. To 
this end we include some possibly helpful identities in the appendix, as well as a discussion of the issue of 
multiple poles. 

We close with a discussion (Section 5) of the tree- level input and possible alternative approaches. We 
compare our approach of spinor integration to the procedure of Ossola, Papadopoulos and Pittau of 
evaluating integrands at particular values of loop momentum and solving a system of linear equations. 
Both our method and theirs algebraic substitution at the integrand level. We claim that the complexity is 
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equivalent. 



2. Unitarity cuts with massive propagators 

Here we develop the program of double-cut phase space integration in (4 - 2e) dimensions (T|, 0, 



19, ^] by generalizing to the case of propagators with nonzero and different masses. 
One difference from the massless case is in the basis of master integrals. If all propagators are massless, 
the basis consists of scalar pentagons, boxes, triangles and bubbles. But if some propagators are massive, 
there are tadpoles and massless bubbles (i.e., the case K 2 = for bubbles). The tadpoles are 16] 



2 

= !^_ +m 2( 1 _ 7 _2iog( m )) + 0( e ), (2.2) 

while for massless bubbles we find 

M*- 2£ -Mt 2£ T[l + e] 
h{K ~ 0) ~ Ml - Mf e(l - e) (2 ' 3) 

= \ + (1 - 7) - M 2 2 _ M 2 ( M l lQ g( M i) " M * log(M 2 )) + 0(e) (2.4) 

or - - (7 + 2 log Mi) + 0(e) if Mi = M 2 . 

It is evident from the lack of all momentum dependence that unitarity methods alone can never detect 
these two functions. 

However, as discussed in [jl6f| , it is possible in some cases to address this difficulty by considering the 
known divergent behavior of the amplitude. For massless gauge theory amplitudes, the conditions are 
such that quadratic divergences cancel and the remaining (logarithmic) divergence matches the known 
value. With these conditions we can fix coefficients at leading order in e. In more general cases, we might 
need some additional information. This could come, for example, from taking the heavy mass limit of one 



propagator [£T|. We shall not discuss the tadpoles and massless bubbles further in this paper. 

2.1 Cuts of Scalar Integrals 

We define the n-point scalar function with non-uniform masses as follows: 1 

r d 4 ~ 2e v 1 
IJMi,M 2 , mi, ... ,m n - 2 ) = -, — * . (2.5) 

This definition is written with a view towards taking the unitarity cut. The masses of the cut propagators 
are Mi and M2, and the momentum flowing through the cut is K. The other momenta Pj and masses rrij 
are the ones necessary to complete the scalar function, and not necessarily in cyclic order. 

1 Note that, for ease of presentation, we are omitting the prefactor i(— l) n+1 (47r) 13 ^ 2 (which was included for example in 
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Our calculations are done in the "Four Dimensional Helicity" (FDH) scheme, i.e. all external momenta 
Ki are 4-dimensional and only the internal momentum p is (4 — 2e)-dimensional. Thus it is useful to write 
p = £ + p,, where i is 4-dimensional and fi is (— 2e)-dimensional. 

d 4-2e p _ f d 4 I f d~ 2e £ e _ f dH (4vr) £ /" d 2/ 2>>-l-e 



(27r) 4 - 2e J (2^) 4 J (27r)- 2e J (2vr) 4 r(-e) 
Then the scalar function as defined in (12. 5D becomes 



I n {Mi,M 2 ,mi,...,m n - 2 )= / dfi 2 ^ 2 )- 1 - 1 



r(-e) J j (2vr)4 

1 



{i 2 _ M 2 _ ^ _ K) 2 _ M 2 _ M 2) U n-2 ({e _ p. )2 _ ^ _ ^ 



(2.6) 



Following the setup of [O], we decompose the 4-momentum into a null component and a component 
proportional to the cut momentum K. 

I=£ + zK, £ 2 = 0, ^\dH=\dz dH5\e)W.K). (2.7) 

We define the "signature of the cut," A[K,M 1 ,M 2 ], as 

A[K, Mi,M 2 ] = (K 2 ) 2 + (Ml) 2 + (Ml) 2 - 2K 2 M\ - 2K 2 M 2 2 - 2M?M$, (2.8) 

and a dimensionless parameter it to be 

4K 2 fi 2 

U - A[K,M u M 2 y (2 ' 9) 



Then 



A[K,Mi,M : ^~ e rl 



2M 1 du u- 1 -". 



4K 2 

Since ( 2 J)4r(-e) ( A ^4j^' M ^ ~) * s a universal factor on both sides (basis and amplitude) of the cut calcu- 
lation, we neglect it henceforth. 

Thus, we have rewritten the integral as 

I n (M 1 ,M 2 ,m 1 ,...,m n _ 2 )= du J dz dH 5 + (I 2 ) 

(21 ■ K)— = 5-^, 

(£2 _ M 2 _ ^ _ jQ2 _ M 2 _ ^2) n n-2 (( ^ _ p.)2 _ m 2 _ ^ 



where /i 2 is related to u through (2.S) 
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2.2 Kinematics and the domain of integration 

In this subsection we determine the integration domain. 

Assuming K 2 ^ 0, we choose the frame where K = (K, 0,0,0), I = (x,y, 0,0). Then the second cut 
propagator is £ — K = (x — K, y, 0, 0). Prom the on-shell conditions, 

x 2 - y 2 = Ml + fi 2 , (x - Kf -y 2 = M 2 + fi 2 , 

we solve for x and y to find 

x = # + ui-H$ i y=± ^_ M? _^ (2 , u) 

The requirement that y has a real solution is the following constraint: 

, 2 <^«1, (2.12) 

with the definition of A[K, M\, M 2 ] given in ( |2.8| ). The condition ( 2.12[ ) also requires the right hand 
side to be positive, which can be arranged by working in a region with sufficiently large cut momentum, 
K > M\ + M2. Then the physical constraint ( |2.12| ) restricts u to lie in a unit interval: 

«6[0,1]. (2.13) 

After using (|2.7| ) to decompose the vector t = £+zK, the lightlike condition I 2 = becomes (x — zK) 2 = y 2 , 
or 

z 2 K 2 - z{K 2 + M\ - Ml) + (Ml + /x 2 ) = 0. 

Solving this equation, we find 



= (K 2 + Ml - Ml) ± y /A[K,M 1 ,M 2 ] - 4K 2 fi 2 
2 2K 2 ' 1 ' ' 

In the kinematic region of this cut, K > 0. Then the positive-light-cone condition 5 + (l 2 ) is equivalent to 
x — zK > 0. Consequently, exactly one of the two solutions ( |2,14| ) for z is selected. Specifically, 

_ (g» + Ml - M 2 ) - y/A[K, Mi, M 2 ] - 4if V 
z- ^2 . ( 2Ab ) 

which we rewrite as 



a-by/l-u 

z = , (2.16) 



where we have defined two useful parameters: 



_ K 2 + M 2 -M 2 \ A A..U;..U, 

a- R2 ,0- K2 ■ {t-Li) 



In the massless limit a = b = 1. 
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2.3 First steps: separating the four-dimensional integral 

Now we are ready to discuss the cut integral. Consider the two delta functions of the cut propagators, 

S(P - Ml - li 2 )8{(£- Kf - Ml - fi 2 ). 

Inside the integral, we perform the following manipulations: 

<$(? - Ml - li 2 )5{(£- K) 2 - Ml - ix 2 ) 

= 5(I 2 - Ml - li 2 )5{K 2 -21-K + Ml - Ml) 

= 5(z 2 K 2 + z(2£ ■ K) - Ml - li 2 )5{{\ - 2z)K 2 — 21-K + Ml — Ml) 

= 5{z{l - z)K 2 + z(Ml - Ml) - Ml - li 2 )5((1 - 2z)K 2 — 21-K + Ml — Ml) 

In the last step we used the second delta function to find 21 ■ K = (1 — 2z)K 2 + Ml — Mf . We will make 
this substitution in the measure ([O]). At this point the first delta function is independent of I. The cut 
part of the integral now takes the form 

J du u^ e J cfe((l - 2z)K 2 + Ml - Ml)5{z{\ - z)K 2 + z(Ml - Ml) - Ml - li 2 ) 

J d 4 £5 + (£ 2 )5((l - 2z)K 2 -2t-K + M%- Ml). 
The second line is the 4-dimensional phase space integration, which can be performed in various ways, as 



discussed in [IS]. We can integrate out z using the delta- function in the first line. Here we need to account 
for the factor 

— 0(1 - z)K 2 + z(Ml - Ml) - Ml - Lt 2 ) = (1 - 2z)K 2 + {Ml - Ml), 



which serves to cancel the first factor of the first line of (|2. 18 ) . Finally, we arrive at the expression 



J du u- 1 - 6 J d 4 £5 + (£ 2 )6((l - 2z)K 2 - 2t ■ K + Ml - Ml), (2.18) 
where z is related to u by ( [2.16 ). 

3. Recursion and reduction formulas with mass 

In this section we compute the massive analogs of the recursion and reduction formulas for master 
integrals presented in IE] and derived in detail in [20|. It is not hard to check that the results of this 



section reproduce the massless results when we set all rrtj = (and hence a = b = 1 in ( 2.17|) ). 

We refer to (2(| rather than review all details of the setup here. However, we do need to remind the 
reader that the results for cuts of these basis integrals were derived with spinor integration, in which the 
massless 4-dimensional vector £ is rewritten as 

I = tXX, (3.1) 
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and the measure transforms as 



<x> 



d A £ S {+) (£ 2 ) (•) = dtt (X,dX)[\,d\](»). (3.2) 
Jo Jx=\ 

Here t ranges over the positive real line, and A, A are homogeneous spinors, also written respectively as 
\£), \£] in many expressions involving spinor products. The first step in spinor integration is to integrate over 
the variable t simply by solving the delta function of the second cut propagator. In this section, we sketch 
the beginnings of certain derivations by writing the integrand before and after this t-integration-which, we 
emphasize, is not true integration. 

Spinor integration proceeds by writing the integrand as a total derivative plus delta functions using 
"holomorphic anomaly" formulas, and finally identifying the contributions of delta functions as residues of 
a complex function. For an exposition of this technique, we refer the reader to || [?J . 

3.1 Bubble 

The cut bubble is exactly the integral described in the previous section without any additional factors. 
We just need to do the four-dimensional integral. 



J dz d 4 £ 5 + (£ 2 )(2£ ■ K)5{z{\ - z)K 2 + z(Mf - Mf ) - Mf - /i 2 )<5((l - 2z)K 2 — 2£-K + M\ — M. 

= J dz d 4 £ S + (£ 2 )((l - 2z)K 2 + Ml - Mf) x 

5(z(l - z)K 2 + z{M\ - Ml) - Ml - n 2 )5((l - 2z)K 2 — 2£ ■ K + Ml — M 2 2 ). 

After phase space integration we are left with 

(1 - 2z)K 2 + Ml - Ml 



K 2 



by/1 - u, (3.3) 



where we have put in z given by (2.16). So, the cut is 



C[I 2 (M 1 ,M 2 ;K)]=b [ du u' 1 '* VT^u, (3.4) 
Jo 



where b is given in ( 2.17| ). This factor of b is where this expression differs from the massless case. 



Recursion relation: Now we define a sequence of functions indexed by nonnegative integers n: 

Bub (n) = C du u- l - e u n ^/T^u~. (3.5) 



o 



Note that the physical cut is related to the zeroth function by a factor of b: 

C[I 2 (M 1 ,M 2 ;K)] = b Bub (0) . (3.6) 

Because the definition ([3l]) is identical to the massless case, we derive exactly the same recursion formula 
as in pi, [p: 

Bub< n > = F^BubW, (3.7) 
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where the form factor is 

(-6)3 



4% = r '\ ■ ( 3 - 8 ) 

(n-ejs 



2 



Here (x) n = T(x + n)/T(x) is the Pochhammer symbol. 

Written in a form suitable for reading master integrands, we have the result 



b du vr l - e u n y/l^u = fQ 2 C[I 2 (M 1 ,M 2 ; K)]. 
Jo 



'0 

3.2 Triangle 

The integrand to start with is 

S(£ 2 - M 2 - fi 2 )5((I- K x ) 2 - M 2 - /x 2 ) 



After t-integration we get 

(, . M? — Mn \ 1 , r- 



(3.9) 



; m^iiiQii] (i\K x \t\(i\Q\e\> 

where 

The four-dimensional integral gives 

1 (1 - 2z) + 1 _ 2 In ^ V ^ , (3.11) 



V^V ^i 2 / V-(2^i-Q)-VA 

where 

A 3 = 4((K 1 .Q) 2 -KfQ 2 ). (3.12) 

The ingredients for ( |3.12| ) are 

K l .Q= Kl + M J 2 - M2 (K l -K,) + {K 2 + M 2 -rnl) 



Q 



2 



(2K 2 (K 2 + M 2 - m 2 ) + (K 2 + M 2 - M 2 )(2K X ■ K 3 )) 2 - (1 - ^A^qA^, M x , M 2 ] 



4(K 2 )3 
where 

A 3;m =o = 4(^1 • K 3 ) 2 - ^i^ 3 2 )> (3-13) 
which we recognize as the signature of a triangle with massless propagators. 
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We can now see that 



A 3 = = ^ (1 - M )A 3 ^ , (3.14) 

It is interesting to see that A3 is built from the factors A[K\, Mi, M%], the signature of the cut, and A3 ;m= o, 
the signature of the triangle. 
Now we define 

— , m^m (3.i5) 

Then we have 

1 , (Z + ^/l-ii 



/•1 1 

C[I 3 (M l ,M 2 ,m 1 ;K 1 ,K 3 )}= / dun' 1 -* (-bVT^H) -= In 



d« „-i- [ -- =J= ) In f Z + . (3.16) 

This is the same expression as in the massless case, but now with a different Z. 
Recursion/reduction relation: We define the integrals 

Tri^(Z) ee f 1 du u-i-^ln ( - + ^E" ) , ( 3 . 17) 
Jo \Z -y/l-Uj 

where the parameter Z is defined as in ( 3,15| ). The physical cut integral is 



C[h(M 1 ,M 2 ,m 1 ;K 1 ,K 3 )] = * TriW(Z). (3.18) 

V^3;m=0 



In this case, our cut triangle functions ( |3.17 ) do depend on the cut-propagator masses M\ and M 2 via Z. 



Apart from the generalized definition of Z, the formula is the same as in the massless case. Therefore, we 
derive the same recursion and reduction formulas as in |20||, namely 



Tri^(Z) = F^%(Z)Tn^(Z) + F 3 ( " ) 2 (Z)Bub (0) , (3.19) 
where the two form factors are given by 

F 3 { %(Z) = ^(l-Z 2 )", (3.20) 



n — e 



s<») (7) _ y 2Z(i - z*y- k r , 21 3 

F ^ {Z) (k-e)i ■ (3 - 21) 



fc=l - 2 

These are functions of the variable Z, defined for a given triangle by ( |3.15| 
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Equation ( |3.19 ) is not yet in the most applicable form. We return to the language of physical cuts by 
including the factor — l/-\f&-z, m =o from ( |3.18 ). The recursion/reduction formula that we need is thus: 



du u~ l - e u n 



1 , fZ + y/l-u 
In 



a/ A 3;m=0 \Z — vl 



u 



F^%(Z)C[I 3 (M 1 ,M 2 ,m 1 ;K 1 ,K 3 )} + F™ 2 (K 1} K Z )C[I 2 {M U M 2 ; A^], 



(n) 



where 



a£ 2 (A 15 A 3 )) 



1 



by A3 ;m= o 
3.3 Box 

The integrand to start with is 

6{I 2 - M 2 - n 2 )5({£- A) 2 - M 2 - fi 2 ) 



((£ - Pi) 2 - m 2 - n 2 )((l - P 2 ) 2 -m 2 - fj, 2 ) 
We define another useful mass-dependent parameter: 

if + M 2 - m 2 



A 2 



After t-integration we get 



b^T 



u 



A 2 {e\Qi\t]{l\Q2\lY 



where 



(3.22) 



(3.23) 



(3.24) 



Qi = {-bVl^i)P l + (at - ^^{a ~ by/l=u)\ A. 

Now the procedure is the same as in the massless case, but with this more general definition of Qi 
Define 

R% = Qi\ u= o 



(3.25) 



-hP, + ( ^ - ^^(« " b) ) A . 



The physical cut is 



C[h(M 1 ,M 2 ,m 1 ,m 2 ;K,P u P 2 ) 



(3.26) 
(3.27) 

(3.28) 



1 



In 



D-Cu + y/T^Il^B - A 



o ' 2A 2 y/B - Au^\D-Cu- sfT^uyjB - Au 
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where 



A 



det 



D 



det 



/ Pi Pi-p 2 pi-k 
Pi ■ P2 Pi P2 ■ k 

\Pi-K P 2 -K K 2 

R 2 R\ ■ R 2 \ 



Rl • R 2 R\ 



K 2 \ Po K K 2 



D = Ri- R 2 . 



(3.29) 



Here again, the form of (|3.28| ) differs from the one in massless case only by the factor of b. 
Recursion/reduction relation: We define 

D - Cu + VT^TlVB - Au 



rW n ) (A, b,c,d) = f 

Jo 



du u 1 6 = In , 

y/B-Au \D-Cu- ^JT^uyfB - Au 



The physical cut is related to the zeroth function in ( [3,30 ) by 

b 



C[h(M 1 ,M2,m 1 ,m 2 ; K, P 1 ,P 2 )] 



2K 2 



Bo^°\A,B,C,D) 



(3.30) 



if A, B, C, D are defined as in ( ^29|) . 



In Table ( 3.31 ) we have listed which kinds of triangles a box with given cut would reduce to (the 
propagator mass must also be correctly identified): 



Box Cut K 


Pi 


P2 


Triangle One's (K\,K 3 ) 


Triangle Two's (K\,K 3 ) 


Ki 


K X2 


-K 4 


(Ki,K 34 ) 


(Ki,K 4 ) 


K 2 


K23 


—K\ 


(K 2 ,K 41 ) 


(K 2 ,Ki) 


K 3 


K 34 


-K 2 


(K 3 ,K 12 ) 


(K 3 ,K 2 ) 


K 4 


K 41 


-K 3 


(K 4 ,K 23 ) 


(K 4 ,K 3 ) 


K12 


K X 


-K 4 


(K 34 ,K 2 ) 


(K 12 ,K 4 ) 


K23 


K 2 


-K x 


(K 41 ,K 3 ) 


(K 23 ,Kx) 



(3.31) 



Boxes are related to triangles and bubbles in the reduction formulas. In these relations we make use of a 
quantity combining variables of the box and the associated triangles to Table ( 3.31| ) : 



C Zi = (Z 2 - l)C + D. 



(3.32) 



Box( n )(^, B, C, D) = F^ 4 {A, £)Box(°)(A, B, C, D) + F^ 3 (A, B, C, D; Z^Tri^Zi) 
+ f[%{A, B, C, D; Z 2 )Tm^{Z 2 ) + fQ 2 (A, B, C, D; ^Bub^ 

where the form factors are given by 

Mi 



F ( 4 %{A,B) 



(n 



£X r, 

~A 



(3.33) 



(3.34) 
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f£ 3 (A5,C,Z);Z,) = ^— X)(fc-l-e)i F£" 3 1) (^). (3-35) 

x^Cfc-l-c^ (f + ^fc^)) • (3-36) 

fe=l V / V 1 2 / 

Again ( 3.33| ) is not the final formula we are after. To get the proper physical result for identifying 
integrands, we need to replace the kinematic factor b/2K 2 . The result is 

1 , -i- t n\ b , (D-Cu + VT^uVB - Au 

du u u . = In ■= 

_2K 2 ^B - Au \D-Cu- VT^uVB - Au 

F 4 ( S 4 (A,B)C[h(M 1 ,M 2 ,m 1 ,m 2 ;K,P 1 ,P 2 )] 

2 

+ £ f[%{A, B, C, D; Zi)C[h{M u M 2 ,mf ; K® , K®)] 
i=i 

+ Fi%(A, B, C, D; Z i )C[I 2 (M u M 2] K)\, (3.37) 

where for the triangles, and K<p are given by Table ( 3.31| ), the mass must also be interpreted 
correctly according to the same table, and the form factors are 



f£! 3 (A,B > C',D;Z < ) = -4ft- Fi%(A,B,C,D;Zi) 



2K 2 

F[%{A, B, C, D; Zi) = F[%{A, B, C, D; Z-). (3.38) 

3.4 Pentagon 

The integrand is 

5Jf - M 2 - y 2 )5((£- Kf - M 2 - /i 2 ) 

{(£- Pi) 2 -m\- [i 2 W- P 2 ) 2 -m 2 - fJ, 2 )((I- P 3 ) 2 - m\ - fi 2 ) ' 

The physical cut is 



where 



/ / \ P 2 + M 2 - mf - 2z(K -Pi) , 

Qi = - (bVT^) Pi + 1 1 ^ -K. (3.39) 



The total integral (apart from the universal prefactor) is 



b / 1 du vT^y/l^u [ {£ dl) [£ M] —^-, j— - jffl^ . m j—= (3.40) 
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Apart from the factor of b and the modified definition of Qi, the analysis proceeds as in the massless 
case p0| . Thus we can cite the result directly as 



C[h{Mi, M 2 , mi, m 2 , m 3 ; if, Pi, P 2 , P 3 )] 

Q 3 • Q2 



.70 



dim 



-1- 



vr 



2^2 



+ 



+ 



S[Qs 


Q2, Qi 


K] 


V(Q 3 ■ 


Q2) 2 - 


Q3Q2 


S[Q 3 


Qi,Q2 


K] 


V(Q 3 • 


Qi) 2 - 


Q3Q1 




Qi,Qs 


,K] 




Q1) 2 - 


Q\Q\ 



In 



Q2) 2 



Q\Q\ 



In 



In 



Q 3 • Q 2 + 

Q 3 • Qi - 7(07 



Q2) 2 



Q3Q2 



Qi 



Q3Q1 



Qi) 2 -QiQl 



Q2-Qi 



Qi) 2 



Q2Q1 



where S[Qs, Q2,Qi, K] is a rational function defined as follows: 



S[Q, 



with 



-8det 



(Qs 
Q2 
\Qt 



K Q 2 -K Qx-K 

Q3 Q\ Q2 ■ Qi 
Q3 Q2 • Qi Q\ 



Qi,Q 3 ,K] = ^, 



-4det 



(Ql 



Q2 ■ Q3 Qi • Q3 
\(h-Q*Q-2-Qi Qi 



Q2 • Q3 Q\ Q2 ■ Qi 



(3.41) 



(3.42) 



(3.43) 



Let us make a few comments on the behavior of pentagon cuts. There are three terms. Each term 
looks like a box signature multiplied by the factor S[»]/(2K 2 ). It is significant that each function S[»] has 
the same denominator I^Hj w hich does not depend on the order of the first three arguments. This can 
be considered as another signature of a cut-pentagon integral. This feature makes the reduction simple. 
Where we see a factor of u n , we just need to write u n S = P{u) + A Y^=\ where P(u) is a polynomial 
in u and A is constant in u. 2 The A term will be the pentagon coefficient, while P{u) indicates reduction 
to boxes. In the massless case, pentagons contribute to terms of 0(e), so they can be neglected. However, 
in cases where propagators are massive or we wish to compute to higher orders in e, their contribution 
must be included. 



4. Formulas for coefficients from double cuts 

Systematic extraction of coefficients of master integrals from four-dimensional spinor integration tech- 
niques has been described in [§,0, building on earlier techniques reviewed in [23]] . When applied to a 
specific amplitude in practice, there are choices to make regarding how to "split" the integrand in partial 
fractions as well as choosing arbitrary spinors on which the final answer does not depend. 

2 It is essential that because all three Si[»\ have same denominator T2, after reduction we have the same A multiplying all 
three of the Si[»]. 
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In this section we present canonical choices to further simplify the method, aiming to automatize the 
heart of this procedure to allow for easy implementation into a computer program. 

For brevity, we do not review the entire spinor integration technique here; instead we refer the reader 
to the explanations given in |], |7]]. The steps that concern us here come after the coordinate change and 
t-integration mentioned at the beginning of Section 3. At this point we have an integrand whose terms are 
rational functions of spinor products and homogeneous in the spinor integration variables. The following 
steps are to split the integrand into partial fractions (using Schouten identities), followed by identification 
of master integrals and integration over a single Feynman parameter. Finally we apply the holomorphic 
anomaly to complete the spinor integration by extracting residues. 

Here we are able to give explicit algebraic functions for coefficients. In the appendix we give a more 
detailed discussion of how to evaluate these functions in practice. 

A summary of the results this section may be found in Section 4.4. 

4.1 Canonical splitting 

Recall that our starting point is the product of the two (on-shell) tree-level amplitudes from each side 
of the unitarity cut. In spinor notation, this integrand takes the general form 3 



Here C is some expression that does not depend on any integration variable, and we have used the relation 
( |2.7| ), i = i + zK, to rewrite the numerator. As we emphasized in the beginning of Section 3, the t- 
integration is trivial because of the second delta function, <5((1 — 2z)K 2 — 2i • K). Thus we can write the 
result of this step immediately: 

(1 ~ 2z)K 2 v / (l-2z)K 2 \ 711 (£\K\trCU-=i (oi ^[^h] 

{i\K\i\ 2 M {t\m ) HSiiWM 

So we see that, after integrating over t, the integrand for phase space integration is a sum of terms of the 
following form: 

term (W 1 nil mil' 

Here G(X) is some monomial in the holomorphic spinor only; hence it factorizes as J| c i) times a constant 
in I. The functions a,j may depend on A as well, for example as [dj\ = (A|Q|. But a,- is the quantity we 
will need in order to split the term further and extract the residues from multiple poles. 
We aim to describe all coefficients in terms of the quantities K, Qi, a>j and G(X). 



3 In this section we deal specifically with the case of massless propgators. The generalization to the massive case is 
straightforward, as in the previous sections. 



_ 14 _ 



The term (|4.l| ) depends as well on n and k. The exponent n is related to the type of master integral 
involved. Terms with n = contribute only to boxes and pentagons. Terms with n = 1 contribute to 
triangles in addition, and terms with n > 2 contribute to all of these plus bubbles. If n = we can multiply 
both numerator and denominator by (^|-ftT|^], so we shall always assume n > 1. 

The first step in our program is to isolate poles by splitting Iterm using the following partial fraction 
spinor identity: 



<£|Qi|<]^|Q2|<l <^|C?2QiK> <^IQi|<] 

Because there are different factors in the denominator, there are different sequences of splitting leading to 
different but equivalent expressions. Here we choose a canonical sequence: 

(1) First we keep untouched and split the factors among themselves. At the end of 

this step each term takes the form - — — 



(e\K\e\"(i\Qi\e\- 

2) Then we split from as often as necessary to get two types of terms, and 



(l\K\£\(£\< 

(3) Finally we will be left with 



(4.3) 



Our task is to find expressions for G/ c (A),i ? j(A) and bj in the most compact and simple form. The 
result is 

m) = ( asaa ^ (4 4 ) 

{(iiKQMr-'nUt^WQtQW ) ' 

C (X) - V G(A) n ^ n "=^" P [aiim v = 2 n US) 

bj = a j+n+k-p- (4.6) 



While -Fj(A) finds its simplest and most compact form in Q4.4D , G P (A) in Q4.5| ) may not be the simplest. 
For example, it can be shown that all terms with 1/ (£|X|£] n can be summed into a single term. In general 
there are several terms with 1/ (£|X|£] a with 2 < a < n. As we discuss in Section 5.3, if we know that 
n < 4, we can work them out explicitly. However, in Section 4.3 we shall use a slightly different method 
to deal with these rational terms. 

The identities necessary for carrying out steps (1) and (2) are the following, which may be proved by 
induction, making use of the basic splitting identity fl4.2|) . 



nti mm U nl Wi w™w 
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n?=iM] n?=i wjg i , w _ P n?=r 2 W) k- p -iI^) ngt> <l 

<^Kr^|QK] (^Ql^" 1 <^K](£|Q|£] ^ l J (£\K\£] P+2 (£\QK\£) n - p - 1 1 ' j 

One necessary condition for the form ((Oj) is that all the Qi are different, which is satisfied for generic 
momenta. 

4.2 Box and triangle coefficients 

In this subsection we derive a closed expression for box coefficients and systematize the derivation of 
triangle coefficients. 

As we know from Section 3.3, box integrals and their coefficients may be labeled by momenta Pi,Pj 
along with the cut momentum K. Prom these momenta we have defined Qi, Qj as in ( p. 25 ). 
The terms we need to consider are of the following form: 

r(i) _ / p(\) 

J l{ '(iiKieMtiQiW 

* dx I {i M) [di df] ( mmm) • R=xQ ^^- 



As discussed in the previous subsection, there will be two types of poles: those from (£\KQ\£) contribute 
to triangles, and those from {£\QjQi\£) contribute to boxes. Let us see precisely how these arise. 
We can construct two massless momenta from Qi and K: 

P[% = Qi + xf 2 K, (4.9) 



where 

I/O TS _L . / A li.\ 

A« = {2Qi ■ Kf - AQ 2 K 2 . (4.10) 



(i) -2Qj ■ K ± \/AB 

X l,2 ~ 2R 2 



(i) 

We make the choice n = P x and find that 



o J \U P 2 W ) [P 2 W pW] (£\R\£] x(x y ;> + 1) - 1 



The poles in this expression come from P«(A) and Y Note as well that since we have chosen n = P® 

the factor [p[ i] £] appears in the numerator and therefore there is no contribution from the pole at 1 1 pjy 



which one might naively expect from a factor of (£\KQi\£) in the denominator of P«(A). 

Apply a partial fraction expansion to the x-dependent factors in the denominator of the integral. The 
result is 

/« - I' dx [ (£ d£) \d£ dA ( m) ( ^ + 1} + i^-m \ \ (4 12) 

~k r } [ £l { mm \ x{x f + i} _ i + x m - m + «\k\£\ ) ) • (4 - 12) 
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where we have used the fact that 



(x\> -xy) 

First term of ( [4.12 ): Let us consider the two terms of ( |4.12|) separately. The first is equal to 

1 r ( wJW ( (J$) i_ -n \ \ 



= f dx [ (£ di) [di d(\ ' FiW 



-i 



o J "\{t\Q,K\ei \ !(!<•■) + 1) _ i_ 

< f<if> [ df *i((ra ( - lnM<>)) 

In the last line we made the substitution (fOj). 

To complete the calculation for I\ , we need to take residue of all poles except the one (£ P x 
from factor (£\KQi\£). This seems difficult to do directly. However, notice that the whole expression is 
holomorphic. Using the result that the sum of residues of all poles of a holomorphic function is zero, we 
get 



\ (i\KQi\i) n rit=wi (*\QtQi\f) 



(4.14) 

These residues contribute to triangle coefficients. In the appendix we show how to evaluate the residue of 



residue of 1 1 



a multiple pole. It can be seen that the expression ( 4.14[) defines the algebraic function for the coefficients 
with input aj and Q. L . 

Second term of ( [4.12 ): The second term of ( |4.12[ ) is defined by 

Again, we do not take the residue of l^£ P®\ from the factor (£\KQi\£). 

There are two kinds of poles in (|4.16| ): one is a simple pole from (£\QjQi\£) and one is a possible 
multiple pole from (l P^/ within (£\KQi\£) n . The former contributes to boxes; the latter contributes 



to triangles. Let us consider the multiple poles first. 

When we replace \£] -» |P 2 (<) ], the spinor dependence cancels out in the integrand, giving 

(£\Qi - K\£] {xf + 1) 

x(£\Q t -K\£] + (£\K\£] ~" x{x f + l) _ x 
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We integrate over the Feynman parameter and find that the residue is 

_( G(X)U^t 2 [^\Qi\i) 



(0 



ln(— x 



2 ) 



(4.17) 



residue of 



3 « v 



Like ( |4.14| ), it will contribute to triangles. 

Now we move to the simple poles in ( |4.16 ). For these simple poles we can evaluate the integral over 
the Feynman parameter first and get 



I. 



(0 



simple poles 



(£ d£) [d£ d £ 



( c(A)n 



tt 2 K\Qi\£) 



In 



(e\K\£\ 



\(i\KQi\i) n T[t=im(t\QtQi\t) WQiWJ ' 



(4.18) 



Now we need to compute the residue from the poles in (£\QjQi\£) ■ Assume without loss of generality that 
i < j, and construct two massless momenta as 



P\ t 2 ~ Qj "^VljiQi) 



(* < j) 



where 



1/1,2 



-2 9 ,.Q^^ A &) = (2Q ( ■ Qj) 2 - iQ]Q]. 



To simplify our expressions, we define the following function: 

G(A)n: = T 2 [a s \Qi\£) 



FiM 



(4.19) 



(4.20) 



(4.21) 



To sum the contributions from these two simple poles (^£ P\ lJ> ^ and {£ P 2 ^y, notice that the P^ 1 differ 
only by a sign in front of the square root. Thus we can expand 



F ( s ) + F M 

i,j T i,j ' 



^ij (-^2 



F : 



(S) 



1,0 



F, 



(A) 



h3 



(4.22) 



Putting this back into the expression fl4.18| ) for the residue, it is straightforward to derive that the contri- 



( (i) 
bution from these two poles in J 2 is 

' ' F S) \u> 



(pPikipP" 


(p 2 ^|Q 8 |P 2 ( ^ 


(p^iq.ipf 


/p? j) \K\P? j) ' 



+ P ( f ) In- 





{p^\K\P^' 


(p[ ij) \Qi\Pi m ~ 


(p 2 ^|Q,|P 2 ^ 



To proceed further, we notice that the same simple-pole factor (£\QjQi\£) shows up in the term 
(or Fj(X) in (4.3)). We can do a similar calculation to get residues in J 2 in terms the function Fji(X). 



One can easily check that P/,i(Pp2^) = Pi,j{P\2 ); by noticing that 



M,2 l-^l"^,! llt=l,t^j,i ^1,2 |ytK2,l 



,(y) 



?(y ) i 



^j(^l 



(4.23) 
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Therefore the sum of contributions from residues associated to (£\QjQi\£) = — (£\QiQj\£) in and I^' is 



Af) 



1 



i 
i 





(p 2 (tf) |0 4 |P a fcr 




(p^lQilPf^ 



+ P<f In 



(p^lOilP^" 


(p 2 (y) |Oi|Pi^ 


(p^IQ^" 


(p^\Q t \P? j) ~ 



F (S) ln 2/t 



(«) . Q 2 



F (A) ln ^ 




'^■(pP) + Pij(P 2 (ii) ) ln 2g, -g f • + p^p^Vp^^) ln Qj/k 2 



-2Qi ■ Qj - VA(y') 



QV K2 



I (.4.26) 



This is the result we are looking for. The first term is the box contribution. To fix sign conventions, recall 
that the double cut is given by 4 



(1 ~ gg) 
K 2 



(1 - 2z) 
K 2 



{£ d£) [i d£] 



1 



(£\Qi\£] {£\Q2\£] 
(l-2z) -1 



tf 2 



-2Qj ■ Qj + 
K 2 -2Qi • Qj - VAW) ' 



:hl 



Therefore the box coefficient is given by 



box;ij 



(1 - 2z) 



(4.27) 



where F^j is given by ( [4.23 ). 

It is important to realize that in fact ( 4.27 ) will be a rational function of u rather than a polynomial, 
because it contains results for pentagons as well. To separate the box coefficient, we need to write 



Cbox;ij{u) = H(u) + AiPi. 

iSpentagons 



(4.28) 



where H(u) is a polynomial in u, Pi is the pentagon cut given by fl3.41|) and A4 is constant in u. Thus 
H(u) will be the true box coefficient, for which we can apply the recursion and reduction formulas, while 
Aj is the true pentagon coefficient. 

The second term in ( [4.26] ) is the final piece needed for triangles. When we rewrite it as 



2 K 2 UeiKQiltril'iilQtQ 



i ln Qi I g(A)n 



n+k-2 
s=l 



a s \Qj 



2 K 2 U£\KQj\i) n ULi 



tQj\ 



residue of (£\QjQi\i) 



residue of {i\QiQj\t) 



4 In this part, we compare to the case where propagators are massless. For the massive case, we need to adjust the definition 
of Qi and include a factors of 6, as explained in section 3. 
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we see that the first term is minus half of the residue contribution of (£\QjQi\£) inside -Fj(A) ( JW) and the 
second term is minus half of the residue contribution of {£\QiQj\£) inside Fj(X) (I^). For the first term, 
when we sum up all simple pole contributions for i^(A), we will be left with the residue of the possible 
multiple pole (£\KQi\£) n . That is to say (using the fact that Qj/K 2 



X i X 



1 ^2 h 



( G(A)n£i- 2 MQ< 



residue of (e\KQ z \i) 



Summing the three triangle contributions ( |4.29| ) with (4.14) and ( 4.17| ) we finally reach 

(-x«) 2 



(4.29) 



1 ( G(X)Uti~ 2 i^U\£) h 

2 \mQi\e) n TlLi,vi(WtQi\e) (W), 

1 ( G(\)ntt 2 iasm£) ln (-4°) 2 N 

2 {mQ^rutx^iiiiQtQ^) {xf x f)j 



residue of 



residue of 



f 5 } 



Notice that 1^4/4) is the signature of triangles. To compare, we recall that for triangles we have 

1 



(1 - 2z) / (£ d£) [£ d£\- 



(£\K\£](£\Qi\£] 
Therefore the triangle coefficient is given by 



(1 - 2z) 



-1 , -2Q i -^ + v / AW 
: m 



\/A» -2Q i -ET-v / AW' 



a 



tri;i 



2(1 - 2z) 



G(X)U 



n+k—2 r i / \ 
8=1 l a s\Qi 



{£\KQ t \£) n Y\ti,t^ (£\QtQ*\£), 

( GWlTZf- 2 [as\Qi\t> 
{(ilKQiWUL^ (£\Q t Qi\£), 



Pi 



(0\ 



(4.30) 



(4.31) 



It is important to observe that both box and triangle coefficients are expressed as the difference of 

(i) (i) 

residues at the two poles and P 2 • This is the origin of the "signature" square roots of box and 
triangle integrals. As we will see shortly, for bubble coefficients, we need to add up the residues at the two 
poles Pi and P 2 . The addition will get rid of the square root signature, but we are left with (1 — 2z), 
which is related to bubbles. 



4.3 Rational part 

Now we address bubble coefficients. They are given by the first term of ( |4.3[ ) with Crfc(A) given by 
fl4.5p . As we have mentioned, for general n and k these are not the simplest and most compact expressions. 
If we constrain n < 4, it is possible to use them directly. Here we use a more general approach. 
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Our starting point is (4.1). We write it here as 

lim 



G(A)nS~ 2 M] 



(4.32) 



where rf = and r is a parameter which we eventually take to zero. 



Because no denominator factor appears more than once, we can use (4.7) directly to reach 

n-1 



E 



G(\m]lt 2 Wj\K + rsr,\£) 



x W K + TS ^ UTA,s^s {*W + rs'ri)(K + rs V )\£) Ui=i {£\Qi(K + rs V )\£) 

i G(\)Uti~ 2 [^\QM) 



+E 



(4.33) 



a ^™ mm uxi (im + T^nM nt ljiVi 

When we take the r — > limit, the second line becomes the contribution to triangles and boxes that was 
discussed in the previous subsection. The first line gives the bubble contribution. It appears that each 
term in the first line diverges in the r —* limit. To proceed, define 

K(s) = K + rsrj. 

Then we find that each term in the first line of ( |4.33| ) can be rewritten as 



1 




G(X) U]^- 2 


aj \K(s)\ij 


U\K{s) - Tsr]\l 


(i\K( s )\e 


T^U\vK{ s )\i) n ~ 2 Yi72 w 


tf-')Ui=i U\QiK(s)\l) 



Now we can expand this expression as a power series in r, keeping those terms that will survive the limit 
r — > 0. Specifically, we substitute 



1 



£\K(s) -tsj]\£ U\K{s)\£ 



h=o (£\K(s) 



Apply the familiar spinor integration procedure to replace one integrand by a total derivative using the 
identity 

[r/ £} n \ ( 1 1 [t] £} n+1 



[£ d£] , 

'"\P\£] n+2 

The result is then the sum of residues of 



(n + 1) (£\P\rj\ (£\P\£] 1 



n+l 



(4.34) 



jn+k-2 „.|^„nA ( n _ 



GWWZi a 3 \K(s)\£ 



n— 1 



T h s h 



h+l 



\h=o (h + l)U\K(s)\£ 



l+h 



(4.35) 



For the multiple pole, (£\rjK(s)\£^ = (£ rj) rj\K(s)\£). Since the numerator contains a factor of [ij £} h+1 , 
we do not take the residue of (£ rj). 
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Finally, after denning the quantity 



R[K(s),{Qi},n]= Res 

poles except i) 



n-l 



r-2 U\ v K(s)\e) UXi^is' - s) UU (l\QiK(s)\£ 



T S 



h=0 {h + l) U\K{s)\t 



i+h 



we can write the bubble coefficients as 



^ n— 1 

C bub bie = —?== y2R[K(s),{Qi},rj\ 



T->0 



The limit is taken by expanding and truncating the series. 

4.4 Summary of results 

Conventions: 

We start from an integrand of the form 



L term 



(i\m n nti (mm 



(4.36) 



(4.37) 



(4.38) 



and present here an expression to be integrated over the final (— 2e) dimensions as discussed in Section 2.3. 
Box coefficients: 



K 2 ( F^P^) + FiM j) ) 



(4.39) 



p[f = Qj + y^Qi (i<j) 



(4.40) 



(ij) _ ~2Qi • Qj ± A (ij)_, „ ,,,2 U - V 2,-Y2 

1 O — "To , i-i 



vi,2 



2Q 



(2Qi ■ QjY - 4QiQ< 



(4.41) 



(e\KQ t \e) n u k t=1 , t ^(i\QtQ 



(4.42) 



Triangle coefficients: 
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Ctri:i 



_ ( G(A)n£i fc - 2 [g«IQ»|l) 



(4.43) 



(4.44) 



(0 

1,2 



-2Qj • # ± \/aW 
2^2 ' 



A» = (2Qi ■ Kf - AQlK 2 



Bubble coefficients: 



(4.45) 



n-l 



Cbubbie = 7= — y~)R[K(s),{Qi},rj\ 
V 1 — M , 



The limit is taken by expanding and truncating the series. 

/ 

R[K(s),{Q i }, V }= Y, Res 

poles except r\ 



n-l 



(4.46) 



/r ^ 2 r h s h (e\ v \i] h+1 



\^=o (/i + l)(^|K(s)|£ 



(4.47) 



K(s) = K + rsr/ 



(4.48) 



For a detailed description of evaluation, see the appendix. 



5. Discussion 

In this section, we first comment on the input needed for this integration program, and then discuss 
the alternative approaches of generalized unitarity cuts and the program of Ossola, Papadopoulos and 
Pittau |]. 
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5.1 Tree level input 

The input needed for the unitary cut method is a collection of (on-shell) treedevel amplitudes, which 
can be calculated by recursion relations. However, results with spurious poles (which generally arise in the 
most compact expressions) can make it hard to apply our method. To avoid this difficulty, we can apply 
spinor identities to regroup the terms into an expression free of spurious poles. 

Let us demonstrate this in one example. The five point function with two massive scalars are given 

by 5 



+ 



(1 2) (2 3) ((h + h) 2 - ^)((l 2 + k 3 )i - M 2)[3|(l + 2)4|1] 
^[1 2f 



K\ 23 \2 3] [31(1 + 2)^11] 

(3|l 2 (l + 2)4|l] 2 

(1 2) (2 3)(1|4|1] <3|£ 2 |3] [31(1 + 2)^11] 
^[1 2]3 



K\ 23 \2 3] [31(1 + 2)^11]" 
The spurious pole in the first term can be split from the others by the identity 

(3|4(1 + 2)|4|1] _ (3|4(1 + 2)|1) (3|4(1 + 2)(1 + 2)|3] 



<1|/ 1 |1][3|(1 + 2)1^11] [3|(1 + 2)|1)(1|4|1] (1|(1 + 2)|3][3|(1 + 2)|4|1]- 
Then 

(3|4(1 + 2)4|1](3|4|2] , (3|4(l + 2)£i|l][12] 



A(£ l ,l + ,2+,3-,£ 2 ) 



+ 



<12)(2 3)[3 2](3|£ 2 |3]<l|£i|l] (1 2) (2 3) [2 3][3|(1 + 2)^1] 
^[1 2f 



K\ 23 \2 3][3|(1 + 2)£ 1 \1] 
Using 

(3|£ 2 (1 + 2)h\l] K\ 23 = (3|£ 2 (1 + 2)^1] [K\ 2 - (3|(1 + 2)|3]) 

= (3|* 2 (1 + 2)€i|l] K\ 2 + (3|(1 + 2)^ |3] [3|(1 + 2)£ 1 |1] + K\ 2 (3|4|1] <3|^ 2 |3] , 

we can add the last two terms together and see that the spurious pole [3|(1 + 2)^i|l] has been canceled (as 
it must). Thus we finally get 

Aft i+ o+ P\ - <3|l 2 (l + 2)li|l] (gjlajgj [1 2](3|(l + 2)l 2 |3) 
[U ' ' ' 2J (12)Kl 3 (l\i 1 \l](3\£ 2 \3] + {12)K1 3 K1 23 " ^ 
A similar calculation gives 

A 9 ^+ P~\ - ^M] , [3 l](2|li|l] (2|l 2 |3] _ (2|l 2 K 123 |2) [1 3] 2 

5il ' ' ' ' 2j (12)^ 3 (3|£ 2 |3](1|4|1] + X^, <3fc|3] K\ 2 K\ 3 K\ 23 ' ^ 

where spurious poles have been canceled. 



5 The sign of the second term is different from the corresponding formula in |24|. We have checked, using the method of 
[e3, that the sign in this earlier formula was a typo. 
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5.2 The quadruple cut 

In four dimensions, the quadruple cut Q is a powerful tool for getting box coefficients because the 
integral is completely localized. Quadruple cuts may be applied here as well since we have separated the 
(— 2e)-dimensional and four-dimensional integrations. We apply them to boxes and pentagons. 

Again, we write in the language of spinor integration (see the brief remarks at the beginning of Section 
3), but this is absolutely not essential to the technique of quadruple cuts. 

Let us start with the quadruple cut (denoted here by C q ) of a box. It is given by 

Cf = J du u- 1 -' J d A l 5 + (£ 2 )5((l - 2z)K 2 -21-K + Ml- M|) 
5((7- P) 2 - M 2 - fi 2 )5((I- P 2 ) 2 - M\ - 

S{Pf + Ml - Ml - z(2K ■ Pi) - 21 ■ Pi)S(Pi + Ml - M\ - z(2K ■ P 2 ) - 21 ■ P 2 ). 
After t-integration, which sets t = -((1 - 2z)K 2 + M\ - Mf)/ (£\K \£], we have 

where 



k\i\ 2 \ (i\K\e] J v {t\K\£\ 



( . s M 2 -M 2 \ Pi + Ml - M 2 - z(2K • Pi) , 
Qi = - ( (1 - 2z) + l — 2 2 j Pi + ^ l - j*f V L K , (5.3) 

To solve the remaining two delta- functions, we need to find two momenta constructed from q^ = Q 2 + XiQi, 
i = 1, 2 such that qf = 0. This is the same construction as in ( |4.20| ). Then the two solutions for the two 
delta-functions are given by 6 

h = \qi)\q 2 ], ^ = |«j)kl]. (5-4) 

Notice that these solutions are complex, as usual for quadruple cuts. There is a universal Jacobian factor, 
which in general is a function of u. Putting it all together we have 

C\ = f du u- l ~ e J{u), (5.5) 
Jo 

where J(u) is the Jacobian and should be symmetric in K,P\,P 2 . 

Next we consider the quadruple cut of pentagons. The calculation is similar, so we shall be brief. After 
t-integration we have 

c« - f 1 du [ a m) u dt\ {l - 2z)K2 + M *- M2 5 5 ( K 2 WQ2\i) \ 

C5 "io J(tdi)[£dl] {im]R2mm ^ {im )6{ {£m ), 



6 These are the of 
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with the Qi defined in ( p.3|) . 

Now we need to sum up the contributions of the two solutions. The Jacobian of the two solutions is 
the same. We end up with 7 

n q f 1 . -!- t J(u)^ f (qi\K\q 2 ] (q 2 \K\qi]\ f 1 -i-e J ( u ) Qfri n n (r. r\ 

Cl= duu 7 |n ; 1 + , , n ; 1 = / duu Tp^2 S [Q2,Qi,Q3,K], (5.6) 

Jo K 2 \m\Q3m\ W2IQ3I91]/ Jo 2K 

where S[Q2, Qi, Q3, K] was defined in (|3l2j ) We see again that the relative factor between boxes and 
pentagons is S[»]/(2K 2 ), just as we learned from the double cut. 

In general, we perform the i-integral first and then add the two solutions. This gives a rational function 
R(u). Then we need to split 

R(u) = P(u) + 2^ a Q3 ^2 ' ( 5 - 7 ) 

Q 3 

where P(u) is a polynomial that gives the box contribution, and oq 3 are constants in u that give corre- 
sponding coefficients of pentagons. This decomposition is the same as the one given by ( [4.28 ). 



After getting box and pentagon coefficients from quadruple cuts, one might continue by applying triple 
cuts to target specific triangle coefficients, and then finally use the usual double cut for the bubble part 



only. Triple cuts have been used to get particular one- loop coefficients in |25|, |26|, |27|, |2q| , and recently a 
nice paper has described a general procedure to compute triple cuts in arbitrary dimensions. One can 
try to systematically study the triple cuts of along the lines presented in this paper. It is easy to see 
that the delta-function there (plus possible derivatives of delta-function) corresponds to our multiple pole. 

5.3 Comparison with OPP method 

Recently, a computation by Ossola, Papadopoulos and Pittau (OPP) ^ has been attracting attention. 
It seems to be a very simple reduction method, which can be performed at the integrand level. The key 
point of the OPP method is that knowing the general form of spurious terms, one can solve algebraically for 
the coefficients of physical and spurious terms from knowing the initial data. The work we have presented 
here is in the same spirit. In our method, it is through splitting into partial fractions that we are able to 
identify contributions to the various basis integrals and find the functions for these coefficients. 

We remark that in it is claimed that in the renormalizable gauge, there is an upper bound on the 
number of spurious terms. For example, there are six for triangles, eight for bubbles. The upper bound in 
the OPP method will correspond to the upper bound of n in our factor (£|iC|^] n . Our experience suggests 
that we will always find n < 4. The reason is the following. Before the t-integration we have 



nf^-^) 2 -^ 2 ) 



7 Here we have inserted the factor 1/2 since we need to sum the contributions of two solutions. This is the same expression 
used in original quadruple cut where on one side, we sum up two solutions and divide by two, while on the other side, we 
have just the Jacobi factor for the basis of box. 
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where the power of i in the numerator N\ has an upper bound, N\ < N 2 + 2. The reason is that in 
renormalizable gauge, the expression derived from Feynman rules has the property that the degree of £ in 
numerator is less than or equal to that in the denominator. After peeling off the two cut propagators, we 
have N\ < N 2 + 2. Now we count the powers of 1/ (£|-fC|£] after t-integration. Numerators contribute —N\ 
while denominators contribute A^. From J t dt and the delta-function, we get an additional —2. In all, we 
have (1/ {£\K\£])~ Nl ~ 2+N2 , thus n = JVi + 2 - N 2 < 4. 

Assuming that indeed n < 4, the correspondence between our method and the OPP method becomes 
clearer. It is easy to see that box and pentagon coefficients have the same expressions. For triangle, we 
have multiple poles up to order 4. From our discussion in the Appendix, it can be seen that a pole of 
order n can be traded for a simple pole with an (n — l)-th derivative. In our case, with n < 4 we have 
up to third derivatives. This is similar to j = 1,2,3 in [^|. For bubble coefficients, we have up to three 
terms, and each one has multiple poles up to order three, so we would expect 3x2 = 6 contributions, 
compared to the eight spurious terms of the OPP method. Our impression is that the spurious terms of 
OPP are hidden as multiple poles as well as simple poles in our language. The algebraic complexity of the 
two methods should be equivalent. 
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A. Evaluation of residues 



In this appendix, we discuss how to evaluate the various functions we defined in section 4 to get 
formulas for coefficients of the basis integrals. Our aim is to make the formulas programmable, so we go 
into some detail. 



A.l Box coefficients 



The box (and pentagon) coefficients are given by (4.27), 



G 



box;ij 



(1 - 2z) 



where 



n+k-2r 



(A.l) 



p(»i)| is\ p(y)" rr A /pfe'lmipfe' 

M,2 l A K2,l LLt=l,t^j,iy 1,2 IVtl-^,! 
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Let us rewrite the coefficient as 



Cbox; 



>.l 



K 2 N 
2(1 - 2z)D' 



where 



N 



n+k-2 k 

N = G{p[ ij) ) H [a. PP](p? j) \K\Pl ij) ] n U (p^\Q t \p[ ij) 

s=l 



n+k-2 

+G(P 2 (ij) ) \~\ [a, p[ ij) ](pi ij) \K\P^ 



s=l 



D (Pl ij) \K\P? j) 



n (^10.1^ 



n (p 2 w) iQ*ipf y) 



For D we use the following identity: 

im\s\m] (V2\s\vi] = {2 V i ■ s)(2 m ■ s) - s 2 (2 Vl ■ m ). 

Then, 



(2K ■ (Qj + y^Qi))(2K ■ (Qj + vT'Qi)) + 
T[Qi,Qj,K]. 



(«')/ 



(A.2) 



(A.3) 



The function T[Qi,Qj, K] in this formula is T2/Q 2 , with the definition of T2 given in ( 3.43Q . Using this 
definition we have 



D = T[Qi,Qj,K] n Yl T[Qi,Qj,Q t ]. (A.4) 

It is easy to see that each term in N may be written as a product of the form 

[a a P2] (PAQ\Px] (Pi b s ) = [a s \P 2 QPi\b s ) = [a s \(Q J +y 2 Qi)Q(Q J +yiQi)\b s ) . 
Thus we see that the evaluation of box coefficients ( |A,2j ) is straightforward. 

A.2 The residue of multiple poles 

For triangles and bubbles, we need to know how to get residue of multiple poles. There are several 
ways to do this. The first one is to use the splitting method discussed in @. 



Our functions \P^ j) ) |P 2 (ii) ] and ^"^P? 3 '] showing up in N and D are the £3, Ia of | 
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The second method is to shift momentum parametrically and then take the limit where the parameter 
goes to zero. We illustrate this method with the example of a double pole. We denote the small parameter 
by r. 



W V? nil (* h) ~* V V) {£ (V + ra)) \[ k (£ bj) 



:r- 

1 U k j=i(vaj) 1 n)=i(( ? ? + ™) a ji 



(V (V + ™)> ]J k ( v bj) (V (V + ra)) ft* (fo + Ta ) b 3 



i nU(na j )( 1 iC=i(i+^: 



■j=n* 1 ' (v b s ). 



T(va)n k j=1 {r,bj) y U k j=A 1 + T 

1 I! jLi fa Oj) / A (« flj) y- (a 6,) 



3± 

fa a > I\U fa bj) [ fa aj) ' (v bj) 



n-=i fa bj) \fri fa a,-> (v bj) 

We see that this is the correct result. Compared to the first method, we have more terms and one extra 
auxiliary spinor r\. However, in this method, the symmetry among aj's and bi's is explicit. 

Now we can use the above idea to get the general expression for multiple poles: 
1 Ui(*<H) , 1 Ui(^i) 



V VT Uj (* bj) n:=o (* fa + sra)) Uj (i bj) 



E 1 1 Hi (fa + rs a) aj 
rrn-l 

so 



' nr=o,^ so <fa + rs o«) fa + sra)) Uj (fa + rs a) bj) 



y- 1 1 Hi (fa + TS o a ) a i) 

^ fa «) n_1 n:=o, s ^ ( s - s °) TI, (fa + TSoa) 6,) 

y. 1 i a fa *)(i+™ofc$; 



i a^oi) /V i rL(i+r So 



r-n— 1 



y. 1 i n,(i+r g0 |f|; 

(r? a)- 1 n, fa b 3 ) YT'l^s - s ) ^(1 + r So |f , 



We see that what we need to do is to expand the expression inside the parentheses as a series in r, keep 
only the terms up to order r ra_1 . Although this expression has one auxiliary spinor rj, the final result does 
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not depend on it. Without writing out intermediate steps, we reach the following result: 



Wi (V (H 



n-l 

E 

s =0 
Ni 



„n-l 



n 



n—l / 
0,s^s \ s 



E 





e n 

. I<ii<i2— <ijv 1 <k+n—2 q=l 

The third method is very similar to the second method. Upon noticing that 

d 1 (n-l)!^)"- 1 



(It 



n-l 



(ir, 



TS) 



(£r, 



TS) 



(A.5) 



(A.6) 



we can rewrite 



1 Eh (* "« 

(£ v) n EL (t h 



il 



d 

n-l 



TS 



i 



i 



dr 



(ir, 



TS) 



(n-l)\(t s^Uj^bj) 



Now we take the simple residue of ■ 
find that the residue is 

Residue 



ts). That is, we set \l) = |ry — ts), take the derivative, and finally 
1 Hi ((V-™) at) 



dr"- 1 \(n-iy.(r, s^Uj ((v ~ rs) 



(A.7) 



It should not be hard to check that (| 
let us check the case of a double pole: 



1 



is the same as ( |A.7| ) if we identify a = s. As a demonstration, 




(I s)(lr,- ts) ]J k (£ bj ) 



residue 




i rij=i (v - ts a j) 

(v-rs s) fg =1 (ri-rsbj) 

This is exactly the same result given by the shifting technique. 

Knowing how to evaluate the residue of multiple poles, we can discuss how to evaluate coefficients for 
triangles and bubbles given by (4.43) and (4.46). We use the third method. For example, for ( [4.43D we 
can shift K — > K — ts. Then using (|A.6| ) we can transfer it to the residue of a simple pole. We can further 
simplify the sum of these two simple poles using the same technique discussed in our previous subsection 
for box coefficients. After that we can take the derivative and the r — > limit by truncating the power 
scries. 
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